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Abstract. Given a reduced affine algebra A over a perfect field K, we present 
parallel algorithms to compute the normalization A of A. Our starting point is 
the algorith m of [Greuel et al. (2010)], which is an i mprovement of de Jong's al- 
gorithm, see [de Jong (1998)1 [Decker et al. (1999)| . First, we propose to strat- 
ify the singular locus Sing(j4) in a way which is compatible with normalization, 
apply a local version of the normalization algorithm at each stratum, and find 
A by putting the local results together. Second, in the case where K = Q 
is the field of rationals, we propose modular versions of the global and local— 
to-global algorithms. We have implemented our algorithms in the computer 
algebra system Singular and compare their performance with that of the al- 
gorithm of [Greuel et al. (2010)[ . In the case where K = Q, we also discuss the 
use of modular computations of Grobner bases, radicals, and primary decom- 
positions. We point out that in most examples, the new algorithms outperform 
the algorithm of [Greuel et al. (2010)[ by far, even if we do not run them in 
parallel. 



1. Introduction 



Normalization is an important concept in commutative algebra, with applica- 
tions in algebraic geometry and singularity theory. We are interested in comput- 
ing the normalization A of a reduced affine if-algebra A, where K is a, perfect 
field. For this, a number of algorithms have been proposed, but not all of them 
are of practical interest (see the historical account in Greuel et al. (2010)| ). 



milestone is de Jong's algorithm, see [de Jong (1998) Decker et al. (1999)| , which 



is based on the normality criterion of [Grauert and Remmert (1971)| , and which 
has been implemented in Singular (see Decker et al. (2012)|), Macau lay2 (see 
"Grayson and Stillman (2010)|), and Magma (see [Bosma et al. (1997)] ). The al- 



gorithm of Greuel et al. (2010)] {GLS normalization algorithm for short), which 



is also based on the Grauert and Remmert criterion, is an improvement of de 
Jong's algorithm. It is implemented in Singular. The algorithm proposed by 
Leonard and Pellikaan (2003)] and Singh and Swanson (2009)] is designed for the 



characteristic p case. It is implemented in SINGULAR and Macaulay2 and works 
well for small p. 

In view of modern multi-core computers, the parallelization of fundamental al- 
gorithms becomes increasingly important. Our objective in this paper is to present 
parallel versions of the GLS normalization algorithm in that we reduce the general 
problem to computational problems which are easier and do not depend on each 
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other. It turns out that in most cases, the new algorithms outperform the GLS 
algorithm by far, even if we do not run them in parallel. 

We start in Section [2] by reviewing the basic ideas of the GLS algorithm. In 
particular, we recall the normality criterion of Grauert and Remmert. In SectionjSj 
we present a local version of the normality criterion which applies to a stratification 
of the singular locus Sing(^) of A. This allows us to find A by a local-to-global ap- 
proach. Section |4] contains a discussion of modular methods for the GLS algorithm 
and its local-to-global version. Timings are presented in Section [Sj 



2. The GLS Normalization Algorithm 



Referring to Greuel et al. (2010)] and Greuel and Pfister (2007)] for details and 



proofs, we sketch the GLS normalization algorithm. We begin with some general 
remarks. For these, A may be any reduced Noetherian ring. 

Definition 2.1. Let ^ be a reduced Noetherian ring. The normalization of A, 
written A, is the integral closure of A in its total ring of fractions Q{A). We call A 
normal ii A = A. 

We write 

Spec(^) = {P C A I P prime ideal} 

for the spectrum of A and V{J) = {P G Spec (A) \ P ^ J} for the vanishing locus 
of an ideal J of A. If P G Spec(yl), then Ap denotes the localization of A at P. 
More generally, if 5 is a multiplicatively closed subset of A and M is an A-module, 
then S~^M denotes the localization of M at S. 

Taking into account that normality is a local property, we call 
N{A) ^ {P £ Spec(^) I Ap is not normal} 
the non-normal locus of A. Furthermore, we write 

Sing(yl) = {P e Spec(A) | Ap is not regular} 
for the singular locus of A. Then N{A) C Sing(yl). 

Remark 2.2. A Noetherian local ring of dimension one is normal if and only if it is 
regular. See de Jong and Pfister (2000)] Theorem 4.4.9]. 



Definition 2.3. Let A be a reduced Noetherian ring. The conductor of A in A is 
the ideal 

Ca = AnnA(A/A) ^ {a e A \ aAC A}. 

Lemma 2.4. Let A be a reduced Noetherian ring. Then N{A) C V{Ca)- Further- 
more, A is module-finite over A if and only if Ca contains a non-zerodivisor of A. 
In this case, N{A) = V{Ca)- 

To state the aforementioned Grauert and Remmert criterion, we need: 

Lemma 2.5. Let A be a reduced Noetherian ring, and let J C A be an ideal 
containing a non-zerodivisor g of A. Then the following hold: 

(1) If (fi € IIom^(J, J), then the fraction '■p{g)/g G A is independent of the 
choice of g, and Lp is multiplication by f{g)/g- 
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(2) There are natural inclusions of rings 

A C RomAi J, J) = - (gJ ■■aJ)<^AC Q(A), a^ipa, (fi^ — , 
9 9 

where (pa ■ J ^ J denotes the multiplication by a ^ A. 

Proposition 2.6 ( [Grauert and Remmert (1971)] ). Let A be a reduced Noethe 
ring, and let J ^ A be an ideal satisfying the following conditions: 

(1) J contains a non-zerodivisor g of A, 

(2) J is a radical ideal, 

(3) V{Ca) C V{J). 

Then A is normal ifj A = Hom^(J, J) via the map which sends a to (pa- 



Definition 2.7. A pair (J, .g) as in Proposition 2.6 is called a test pair for A, and 
J is called a test ideal for A. 



By Lemma 2.4 test pairs exist iff A is module-finite over A. Given such a 
pair {J,g), the idea of finding A is to successively enlarge A until the normality 
criterion allows us to stop (since A is Noetherian, this will eventually happen in the 
module-finite case). Starting from Ao ~ A, we get a chain of extensions of reduced 
Noetherian rings 

^ = Ao c • • • c c A,c ■■■ c A,n^A. 

Here, ^i+i = Hom^.(Ji, Ji) = ]j{9Ji '-Ai Ji)^ where Ji is the radical of the extended 
ideal JAi, for i>\. Note that {Ji,g) is indeed a test pair for A^: 

Remark 2.8 ( [Greuel et al. (2010)| , Prop. 3.2). Let A be a reduced Noetherian ring 
such that A is module-finite over A, and let A C A' C A be an intermediate 
ring. Clearly, every non-zerodivisor g € A of A is a non-zerodivisor of Q{A). In 
particular, it is a non-zerodivisor of A' . Furthermore, if Ca' is the conductor of A' 
in W = A, then Ca' 2 Ca- It follows that every prime ideal Q € N{A') = V{Ca') 
contracts to a prime ideal P = Q O A G ^{^) = V{Ca)- Hence, if (J, (?) is a test 
pair for A, then P 2 J, which implies that Q 3 V J A' —: J'. We conclude that 
{J',g) is a test pair for A'. 

Explicit computations rely on explicit representations of the Ai as A-algebras. 
These will be obtained as an application of Lemma 2.9 below. To formulate the 
lemma, we use the following notation. Let J C A be an ideal containing a non- 
zerodivisor g of A, and let A-module generators uq = g,ui, . . . ,Us for gJ -.a J be 
given. Choose variables Ti, . . . , Tg, and consider the epimorphism 

1 u 

<^>:A[Tu...,Ts]-^-{gJ :a J), T,^^. 

9 9 

The kernel of $ describes the A-algebra relations on the Ui/g. We single out two 
types of relations: 

• Each A-module syzygy 

aoUQ + aiUi + . . . + asUs = 0, G A, 

gives an element + aiTi + . . . + a^Tg G Ker<I>, which we call a linear 
relation. 

• Developing each product y y , 1 < i < J < s, as a sum y ^ = J^k f^ijk^i 
we get elements TiTj~J2k PijkTk in Ker $, which we call quadratic relations. 
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It is easy to see that these hnear and quadratic relations already generate Ker$. 
We thus have: 

Lemma 2.9. Let A be a reduced Noetherian ring, and let J Q A he an ideal contain- 
ing a non-zerodivisor g of A. Then, given A-module generators Uq = g,Ui, . . . ,Us 
for gj :a J, we have an isomorphism of A-algebras 

1 u- 
A[Tu...,Ts]/R= -{gJ -.a J), T,^^, 
9 9 

where R is the ideal generated by the linear and quadratic relations described above. 



The following result from Greuel et al. (2010)] will allow us to find the normal- 



ization in a way such that all calculations except the computation of the radicals 
\/~Ti can be carried through in the original ring A: 

Theorem 2.10. Let A be a reduced Noetherian ring, let J ^ A be an ideal con- 
taining a non-zerodivisor g of A, let A (1 A' (- Q,{A) he an intermediate ring such 
that A' is module-finite over A, and let J' = V JA'. Let U and H be ideals of A 
and d ^ A such that A' — and J' — \H , respectively. Then 

igf -.A' J') = J {dgH -.A H) C Q(A). 

Remark 2.11. In the case where A = K[Xi . . . , is a reduced affine algebra 

over a field K, let Pi, ... , Pr be the associated primes of the radical ideal /. Then 



A-K[X,,...,X„]/P, X ... xX[Xi,...,X„]/P,, 
and A is module-finite over A by Emmy Noether's finiteness theorem (see 



bwanson and Huneke (2006)| ). Thus, using techniques for primary decomposition 



as m Greuel et al. (2010) Remark 4.6], the computation of normalization can be 
reduced to the case where A is an affine domain (that is, / is a prime ideal). When 
writing our algorithms in pseudocode, we will always start from a domain A. Talk- 
ing about a non-zerodivisor then just means to talk about a non-zero element. 

Remark 2.12. If A is an affine domain over a perfect field K, we can apply the 
Jacobian criterion (see Eisenbud (1995)] ): If M is t he Jacobian ideap] of A, then 



M is non-zero and contained in the conductor Ca (see Greuel et al. (2010) Lemma 
4.1]). Hence, we may choose -y/M together with any non-zero element g of a/M as 
an initial test pair. Implementing all this, the GLS normalization algorithm will 
find an ideal J7 C A and a denominator d £ Ca such that 

A=]ucq{A). 

a 

Since Al is contained in Ca, any non-zero element of M is valid as a denominator: 
If 7^ c e M, then c ■ U' is an ideal of A, so that \U = \U' . 

For the purpose of comparison with the local approach of the next section, we 
illustrate the GLS algorithm by an example: 



^The Jacobian ideal of A is generated by the images of the c X c minors of the Jacobian matrix 
( J^), where c is the codimension, and /i, . . . , /r are polynomial generators for /. 
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Example 2.13. For 

A = K[x, y] = K[X, Y]I{X^ + Y^Y - if), 
the radical of the Jacobian ideal is 

and we can take g a; G J as a non-zerodivisor of A. In its first step, starting 
with the initial test pair {J,x), the normalization algorithm produces the following 
data: 

{/« -.^xJ -.A J ^{x,y{y^lf)^ and A^ := A[h] -.^ A[T^]/ h = 
with relations and isomorphism given by 

h = {-Tix + y{y - l)^ riy(y - 1) + x\tI + x^{v - 

and 

^ y{y-lf 
ti i-> , 

X 

respectively. In the next step we find 



Ji \J{x,y{y- I)) A, = (2;, - ^)M)a, 

= ^{x\xy{y-l),y{y~lf)^^: -H,. 

X ^ X 



Using the test pair (Ji, x) and applying Theorem 2.10 and Lemma 2.9 we get 

-(xJi -.A, Jl) = ^{x^H, -.A ffl) 



a; " a;^ 



\{x',xyiy-l),yiy-ir)^^:\u^'^ 



and 



with relations and isomorphism given by 
I2 = {T2X ~ n{y - 1),-T3X + y{y - 1), T2y{y - 1) + x^, T2y\y - if + T^x^ 
Ti + {y^l),T2T3 + x.T^-T2y) 

and 

y{y-lf y{y-l) 
a;^ X 

respectively. In the final step, we find that A2 is normal, so that A = A2. 

3. Normalization via Localization 

In this section, we discuss a local-to-global approach for computing normaliza- 
tion. Our starting point is the following result: 

Proposition 3.1. Let A be a reduced Noetherian ring. Suppose that the singular 
locus Sing(j4) = {Pi, . . . , Ps} is finite. For i = 1, . . . , s, let Si = A \ Pi, and let an 

intermediate ring A C v4('') cAbe given such that = S^'^A. Then 
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Proof. We will show a more general result in Proposition |3.2| below. □ 

That Sing(A) is finite is, for example, true if A is the coordinate ring of a curve. 
Whenever Sing(yl) = {Pi, . . . , Pg} is finite, the proposition allows us to find A 
by normalizing locally at each Pi using Proposition |3.3| below, and putting the 
local results together. In the case where Sing(yl) is not finite, working just with the 
(finitely many) minimal primes in Sing(A) will not give the correct result. However, 
it is still possible to obtain A as a finite sum of local contributions: The idea is 
to stratify Sing(A) in a way which is compatible with normalization. For this, if 
P e Sing(A), set 

Lp= n p- 

PDPeSing(A) 

We stratify Sing(A) according to the values of the function P i-> Lp. That is, if 

C = {Lp\Pe Sing(A)} 

denotes the set of all possible values, then the strata are the sets 

VL = {Pe Sing(A) \Lp^L},Le C. 

We write Strata(^) = {Vl \ L Q £} for the set of all strata. If Pi, Pr denote the 
minimal primes in Sing (A), we have 

/:c{aer^. I rc {!,..., r}}. 

Hence, the set of strata is finite. By construction, the singular locus is the disjoint 
union of all strata. For V E Strata(yl), write Ly for the constant value of P i-> Lp 
on V. 



We can now state and prove a result which is more general than Proposition 3.1 



Proposition 3.2. Let A be a reduced Noetherian ring with stratification of the 
singular locus Strata(A) — {Vi, Vs}. For i — 1, . . . , s, let an intermediate ring 
A C CA be given such that S-^A^^'> = S-^A for each S ^ A\P, P £ Vi. 
Then 

1=1 

Proof. By construction, B := X]i=i ^^^^ ^ ^- We wish to show equality. It suffices 
to show that if P e Spec(^) is a prime ideal and S = A\P, then S~^B = S^^A. 
If P e Sing(A), then P e V, for some i. Hence, S'^A'^^^ = S-^A, and the local 
equality is obtained from the chain of inclusions 

5'-i^(i) c S-^B c S-^A = 

li P ^ Sing(^), then S^^A is normal, and the local equality follows likewise from 
the chain of inclusions 



S^'A c S^'B c S-'A ^ S-^A. 

□ 

For a given stratum V = Vi, the modification of the Grauert and Remmert 
criterion below will allow us to find a ring A^*-' as above along the lines of the 
previous section: 



PARALLEL ALGORITHMS FOR NORMALIZATION 



7 



Proposition 3.3. Let A he a reduced Noetherian ring such that A is module-finite 
over A, and let A Q A' Q A he an intermediate ring. Let V £ Strata(yl), and let 
J' = \J Ly A! . Suppose that Ly contains a non-zerodivisor g of A. If 

A' '^}iomA'iJ',J') 

via the map which sends a' to ipa' , then the localization S^^A' with S = A \ P is 
normal for each P (z V . 



Proof. The assumption and [Eisenbud (1995) Proposition 2.10] give 



S-^A' = S-\}iomA'{J',J'))^lloms-iA'iS'^J',S-^J'). 
Hence, the result will follow from the Grauert and Remmert criterion (Proposition 



2.6[ ) applied to S~^A' once we show that the localized ideal S~^J' satisfies the three 
conditions of the criterion. First, since forming radicals commutes with localization, 
S^^J' is a radical ideal. Second, the image of g in S~^A' is a non-zerodivisor 
of S-'^A' contained in S'^J'. Third, we show that V{Cs-ia') N{S-^A') C 
V{S-^J'). For this, we first note that 

V{Cs-ia) - N{S-^A) = {S-ip I P e N{A), P c P}. 

Indeed, prime ideals of S^^A correspond to prime id eals of A contained in P. 
Let now Q e N{S^^A'). Then, as shown in Remark 2.8 Q contracts to some 



S-^P C S-^A with P e N{A), P <ZP. This implies that 



Q 2 y{s-^P)is-^A') = ^s-HPA') = s-\ypA') d s-V, 

as desired. □ 
In the situation of Proposition [3^ let a non-zerodivisor g € Ly of A be known. 



Then, using {Ly, g) instead of a test pair as in Definition 2.7 and proceeding as in 
the previous section, we get a chain of rings 

A c Ai c • • • c v4„ c 3 

such that 5'~^(j4„j) is normal and, hence, equal to S^^A = S^^A for all S* ^\ P, 
P€V. 

Summing up, we are lead to Algorithms [T] and [2] below. 



Algorithm 1 Normalizing the localizations 

Input: An affine domain A = K[Xi, . . . , A"„]// over a perfect field K, a stratum 

V e Strata(A), and Oy^geLy. 
Output: An ideal U C A and d e A with Q A and {^U) ^ S-'^A for all 

S = A\P, P eV. 

1: return the result of the GLS normalization algorithm applied to {Ly,g); 
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Algorithm 2 Normalization via localization 

Input: An affine domain A = K[Xi, . . . , Xn]/I over a perfect field K. 
Output: An ideal U C A and d e A such that A^ C Q{A). 

1: J := Vm, where M is the Jacobian ideal of A; 

2: choose ^ g & J; 

3: compute the strata of the singular locus Strata(yl) ~ {Vi, ...,Vs}; 
4: for all i do 

5; apply Algorithm [T] to {Vi,g) to find an ideal Ui ^ A and a power di — g' 

with A C j-Ui C A and S^^{j-U,) = S-^A for all S ^ A\P, P e V^; 
6: m := max{mi, . . . ,toJ, d := g"\ [/ g™"™'{7,; 
7: return {U,d); 



Remark 3.4. In Algorithm [2j it may be more efficient to choose possibly different 
non-zero elements gi G Ly.. In Step 5, the algorithm computes, then, pairs {U-,di) 
with ideals C// C A and powers di = g^\ As explained in Greuel et al. (20I0)[ 
Remark 4.3], starting from the {Ul,di), we may always find a denominator d £ M 
and ideals Ui C A such that 2^i — -^U^ for all i. Then, the desired result is 

{T.^U^.d). 

Remark 3.5. In Algorithm[2j it is sufficient to consider the minimal strata, that is, 
the strata V such that Ly is minimal with respect to inclusion. Denote, as above, 
the minimal primes of the singular locus of A by Pi,...,Pr. We can obtain the 
minimal Ly as all possible intersections Plier-^*' with subsets F C {I, ...,r} which 
are maximal with the property that X]igr^« ^ 

Example 3.6. We come back to the coordinate ring A of the curve C with defining 
polynomial f{X,Y) = X*+Y'^{Y-lf from Example [2I3 



to discuss normalization 



via localization. The curve C has a double point of type At, at (0, 0) and a triple 
point of type Eq at (0, 1). We illustrate Algorithm[2] using for both singular points 
the non-zerodi visor g — x: For the ^3 -singularity, consider 

Fi = {x,y)^ and Si=A\Pi. 



The local normalization algorithm yields S-^ ^ A = ^(j^Ui), where 

di = x^ and Ui = (a;^ y{y - I)^)^ . 
For the i?6 -singularity, considering 

P2 = {x, y - 1)^ and S2 = A\ F2, 



we get A = 5^^(^f/2), where 



d2 = 



and U2 = (x\ xy" (y - 1) , y" {y - I)')^ 



Combining the local contributions, we get 



\u = ^U^ + ]-U2. 
d di d2 



U = (x\ xy" (2/ - 1) , y(y - 1)^, y^ (y ~ I)^^^ 
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A moment's thought shows that U coincides with the ideal t/'^-* found in Example 

The local-to-global approach is usually much faster than the global algorithm 
even when not run in parallel. The reason is that the minimal primes of the singular 
locus are much simpler than the singular locus itself. Therefore, in the local-to- 
global case, the intermediate rings are much easier to handle. Most notably, the 
representations of the intermediate rings as affine rings involve considerably less 
variables than in the global case. In the following example, we exemplify this 
difference. 

Example 3.7. Consider the projective plane curve defined by the polynomial 

/i,4 = [x^ + + z'^f - A{x^Y^ + a:^z^ + Y^z^) e q[a:, y, z], 

which will be reconsidered in Section[5]with respect to timings. After the coordinate 
transformation Z H> iX ~ 2Y + Z, all singularities of the projective curve lie in the 
affine chart Z ^ 0. Write 

/ = /i,4(X, r, 3X-2Y + l)e Q[X, Y] W 

for the defining polynomial of the affine curve, and let A = Q[x,y] = W/{f). 

The curve has 15 singular points: the radical of the Jacobian ideal M decomposes 

as 

y/M = {y, 121x^ + U2x^ + 64:X^ + 13x + 1) n (y, 2x + 1) 
n {211y'^ - 131y^ + 51y^ - Uy + l,3x-2y+ 1) 
n (lly^ - 23y^ + 1 V - 7y + I, x) n {y + I, x + 1) n {3y ~ 1, x). 

more complicated (see below) . We compare the global approach to the local strategy 
at the singularity corresponding to the test ideal J = (y, 2x + 1). 

In the local setting, we use the non-zerodivisor g — y and compute the ideal 
quotient 

= (y, 292822:^ + 83369x® + 105668x'' + 78296x6 + 37382x^ + 11926x'' 
+ 2542x^ + 349x2 + 28x + 1). 

We observe that in addition to y, the ideal Ui requires only one more generator. 
Hence, the representation of = g (5"^ • J) affine ring Ai = A[Ti]/Ii requires 
only one additional variable Ti . The ideal Ii is generated by 10 linear relations and 
one quadratic relation. Next, we compute the radical of the image of J in Ai. 
Technically, this means to compute the radical \/ J + Ii in the polynomial ring 
M^[Ti] (here, by abuse of notation, we denote the preimage of J in the polynomial 
ring also by J). The ideal Ii is quite complicated. Since J is generated by linear 
polynomials, however, the reduced Grobner basis of J + /i is very simple. It is 
easily computed as 

J + Ii = (y, 2X + 1, 16384Tf -Ti+ 625). 
As a consequence, the computation of the radical 

Ji = y/J + h = {Y, 2X + 1, 128Ti - 25) 
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is cheap as well. In the next step, A2 can be represented as an afhne algebra over 
A with, again, only one new variable. Hence, verifying that A2 is already a local 
contribution to A at the singularity corresponding to J is also cheap. 

In contrast, the global approach uses the test ideal J — \fM, which is generated 
by one polynomial of degree 3 and three polynomials of degree 6. As a non- 
zerodivisor, we consider the lowest degree generator g = 3x^y — 2xy^ + xy. As 
in the local case, the first ideal quotient gJ : J is easily obtained. However, in 
addition to g, it requires three more generators. Hence, as an affine algebra over 
A, it is represented as Ai = A[ri, T2, TaJ/Ji, where the ideal of relations Ii C 
A[Ti,T2,T3] is generated by 6 linear and 6 quadratic relations. No significant 
reduction occurs in J + Ii since J does not contain any linear polynomial. The 
complexity of Buchberger's algorithm grows doubly-exponentially in the number 
of variables. Compared to the local case, this increase in complexity makes the 
computation of ^/J + /i considerably more expensive. In fact, Singular does not 
compute the corresponding Grobner basis within 2000 seconds. 



4. Modular Methods 

Algorithm [2] from Section [3] is parallel in nature since the computations of the 
local normalizations do not depend on each other. In this section, we describe a 
modular way of parallelizing both the GLS normalization algorithm from Section 
[2] and the local-to-global algorithm from Section |3] in the case where K — Q is 
the field of rationals. One possible approach is to just replace all involved Grobner 
basis respectively radical computations by their modular variants as introduced by 
Arnold (2003)| and Idrees et al. (2011)] . These variants are either probabilistic 



or require rather expensive tests to verify the results at the end. In order to re- 
duce the number of verification tests, we provide a direct modularization for the 
normalization algorithms. The approach we propose requires, in principle, only 
one verification at the end. In the local-to-global setup, however, it is reasonable 
to additionally handle the Grobner basis computation for the Jacobian ideal, the 
subsequent primary decomposition, and the recombination of the local results by 
modular techniques. We exemplarily describe the modularization of the GLS nor- 
malization algorithm as outlined in Section [2] Each of the local normalizations in 
Algorithm [2] from Section [3] can be modularized similarly. 

Fix a global monomial ordering > on the semigroup of monomials in the set of 
variables X — {Xi, . . . , X„}. Consider the polynomial rings W = Q[X] and, given 
an integer N>2,Wn = {Z/NZ)[X]. If T C or T C is a set of polynomials, 
then denote by LM(T) {LM(/) \ f £ T} its set of leading monomials. If I e Q 
withgcd(a,6) = 1 andgcd(6,7V) = 1, set (f)^ := {a + NZ){b + NZ)'^ e Z/NZ. If 
f Cz W is a polynomial such that N is coprime to any denominator of a coefficient 
of /, then G Wn is the polynomial obtained by reducing each coefficient modulo 
N as just described. If iJ = {hi, . . . ,ht} is a Grobner basis in W such that N is 
coprime to any denominator in any hi, then write Hj\j ~ {{hi)pf, . . . , {ht)[^}. Given 
an ideal / C W^, set /at = (/at | / e / n Z[X]) C Wn and {W/I)n = Wn/In- 

Remark 4.1. For practical purposes, the ideal I <ZW is given by a set of generators 
fi, . . . , fr- Then for all but finitely many primes p, the ideal Ip can be realized via 
the equality 

Ip^{{fl)p,-Afr)p)^Wp. 
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When performing the modular algorithm described below, we reject a prime p if 
the ideal on the right hand side is not well-defined. Otherwise, we work with this 
ideal instead of Ip. The finitely many primes where the two ideals differ will not 
influence the result if we apply error tolerant rational reconstruction (see Remark 



4.2). 



From a practical point of view, we work with ideals of the polynomial ring W 
containing /, but think of them as ideals of the quotient ring A = W/I. Therefore, 
we simplify our notation as follows: If I C J C W are ideals, then we denote the 
ideal induced by J in A also by J. Vice versa, if J C A is an ideal, then its preimage 
in W is also denoted by J. Similarly for Wn- 

From now on, / = (/i, . . . , fr) C W will be a prime ideal. We wish to compute 
the normalization of the affine domain A = W/I using modular methods. For this, 
we fix a polynomial d G W which represents a non-zero element in the Jacobian 
ideal M of A. This element of M will also be denoted by d. It will serve as a 
"universal denominator" for all normalizations in positive characteristic as well as 



for the final normalization in characteristic zero (see Remark 2.12 for the choice 
of denominators). In characteristic zero, we write U{0) for the ideal of A which 
satisfies 2^/(0) = A, and G(0) for the reduced Grobner basi^of U{0). Furthermore, 
we write V{Q) C A[Ti, . . . ,Ts] for the ideaj^of relations on the elements of ^G(O) 



which represents A as an ^-algebra as in Lemma 2.9 We denote the reduced 
Grobner basis of V{0) by R{0). In the same way, if p is a prime number which does 
not divide any denominator in the reduced Grobner basis of / and such that Ap is 
a domain and dp is non-zero and contained in the conductoiQof Ap, we use U{p), 
G{p), V{p), and R{p) in characteristic p. 

Note that G{0)p is not necessarily equal to G{p). However, as we will show in 
Lemma |4. 5| below, equality holds for all but finitely many primes p. Relying on this 
fact, the basic idea of the modular normalization algorithm can be described as 
follows. First, compute the Jacobian ideal M of A and choose a polynomial d G W 
representing a non-zero element d € M. Second, choose a set V of prime numbers 
at random, and compute, for each p £ V, reduced Grobner bases G{p) C Wp such 
that -J- {G{p)) C Q{Ap) is the normalization of Ap. Third, lift the modular Grobner 
bases to a set of polynomials G C W and define U := (G). We then expect that 
U = U{0) and G = G(0). 

The lifting process has two steps. First, assuming that all LM(G(p)), p G V, 
are equal, we can lift the Grobner bases in the set GV '■= {G{p) \ p G V} to a 
set of polynomials G{N) C Wn, with N := Wp^-pP- For this, apply the Chinese 
remainder algorithm to the coefficients of the corresponding polynomials occurring 
in the G(p), p G V. Second, compute a set of polynomials G C ly by lifting 
the modular coefficients occurring in G{N) to rational coefficients as described in 
Bohm et al. (2012)| : 



Recall that reduced Grobner bases are uniquely determined. For practical purposes, however, 
we do not need to reduce the Grobner bases involved since the lifting process described below only 
requires that the result is uniquely determined by the algorithm. 

''with respect to ideals of W[Ti , . . . , Ts] and A\T\ , . . . , Ts] , we use the same setup and notation 
as for ideals of W and A. 

^From a practical point of view, we check whether dp is in the Jacobian ideal of Ap. 



PARALLEL ALGORITHMS FOR NORMALIZATION 



12 



Remark 4.2. Rational reconstruction via the Chinese remainder theorem and Gauss- 
ian reduction is error-tolerant in the following sense: Let Ni and N2 be inte- 
gers with gcd(A^i,iV2) = 1, and let f € Q with gcd(a, 6) = gcd(A^i,&) = 1. Set 
ri := (1)^ S Z/NiZ, let r2 € 'E/N2I' be arbitrary, and denote by r the image of 
(ri,r2) under the isomorphism 

Z/iViZ X Z/N2Z Z/{NiN2)Z. 

Lifting r to a rational number via Gaussian reduction will generate, starting from 
(ao,&o) = {NiN2,0) and (ai,fei) — (r, 1), a sequence of rational numbers {ai,bi) 
obtained by setting 

{ai-2,bi-2) = qi{ai-i,bi^i) + {ai,bi), 

where qi is chosen such that {ai,bi) has minimal Euclidean length. Computing 
this sequence until the Euclidean length does not decrease strictly any more, we 
obtain a tuple {ai,bi) with ^ = |, provided that N2 <^ Ni. For details, see 
Bohm et al. (2012)| . 



Just as for QV, we proceed for the set of reduced Grobner bases TZV := {R{p) \ 
p G V} giving the modular algebra relations. 

As for other modular algorithms based on Chinese remaindering, we need suit- 
ably adapted notions of a lucky prime and a sufficiently large set of lucky primes: 

Definition 4.3. Using the notation introduced above, we define: 

(1) A prime number p is called lucky for A if U{0)p — U{p), V{0)p — V{p), 
and the following hold: 

(a) Ap is a domain. 

(b) dp is a non-zero element in the conductor of Ap. 

(c) LM(G(0)) =LM(G'(]5)). 

(d) LM(i?(0)) = LM(i?(p)). 
Otherwise p is called unlucky for A. 

(2) A finite set V of lucky primes for A is called sufficiently large for A if 



p > max < 2 



\c\' 



c a denominator or numerator of a 
coefficient occurring in G{0) or i?(0) 



Remark 4.4. A modular algorithm for the basic task of computing Grobner bases is 
presented in [Arnold (2003)| and Idrees et al. (2011)] . In contrast to our situation 



here, where we wish to find the ideal U{0) by computing its reduced Grobner 
basis G{0), Arnold's algorithm starts from an ideal which is already given. If p 
is a prime number, J C is an ideal, H{0) is the reduced Grobner basis of J, 
and H{p) the reduced Grobner basis of Jp, then p is lucky for J in the sense of 
Arnold if LM(7J(0)) = LM{H{p)). It is shown in [Arnold (2003)] Thm. 5.12 and 
6.2] that if p is lucky for J in this sense, then H{0)p is well-defined and equal to 
H{p). By [Arnold (2003) Cor. 5.4 and Thm. 5.13], all but finitely many primes are 



Arnold-lucky for J. Moreover, if P is a set of primes satisfying Arnold's condition 
LM(i?(0)) = lM{H{p)) for all p G and such that V is sufficiently large with 
respect to the coefficients occurring in H{Q), then the H{p), p E V, lift to H{0). 

In our situation, if p is a prime number, we find U{p) on our way, but U{0)p is only 
known to us after U{0) has been computed. Therefore, the condition U{0)p — U{p) 
in our definition of lucky can only be checked a posteriori. Similarly for V{0)p — 
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V(p). However, when performing our modular algorithm, by part 1 of Lemma 4.5 
below and Remark |4.2[ there are only finitely many primes not satisfying these 
conditions and these primes will not influence the result of the algorithm. 

Lemma 4.5. With notation as above, we have: 

(1) All but a finite number of primes are lucky for A. 

(2) If V is a sufficiently large set of lucky primes for A, then the reduced 
Grobner bases G{p), p € V, lift to the reduced Grobner basis G(0). In 
the same way, the R{p), p €V, lift to i?(0). 

Proof. With respect to part 1, it is clear that conditions (la) and (lb) in our 
definition of lucky are true for all but finitely many primes. Moreover, ■^U{0)p 
is integral over Ap for all but finitely many p. Since testing normality via the 
Grauert and Remmert criterion amounts to a Grobner basis computation, and since 
reducing a Grobner basis modulo a sufficiently general prime p gives a Grobner 
basis of the reduced ideal, we conclude that U (0)p = U{p) for all but finitely many 
primes. Furthermore, if U{0)p — U{p), condition (Ic) from our definition of lucky is 
equivalent to asking that p is lucky for U{0) in the sense of Arnold, so that also this 
condition holds for all but finitely many primes. For V{0)p — V{p) and condition 
(Id), we may argue similarly since finding the ideal of algebra relations amounts to 
another Grobner basis computation. 

For part 2, let 7^ be a sufficiently large set of lucky primes for A. Then, as pointed 
out above, G(0)p is well-defined and equal to G{p) for aWp g V. Furthermore, since 
V is sufficiently large, the G{0)p, p G V, lift to G(0). In the same way, we may 
argue for the relations. □ 

From a theoretical point of view, the idea of the algorithm is now as follows: 
Consider a sufficiently large set V of lucky primes for A, compute the reduced 
Grobner bases G{p), p £ V, and lift the results to the reduced Grobner basis G(0) 
as described above. 

From a practical point of view, we face the problem that the conditions (Ic), 



(Id), and (2) from Definition 4.3 cannot be tested a priori. To remedy the situation, 
we proceed in a randomized way. First, we fix an integer t > 1 and choose a set 
of t primes V at random. Second, we delete all primes p from V which do not 
satisfy conditions (la) and (lb). Third, we compute GV — {G{p) \ p G V} and 
TZV = {R{p) \ p &V}, and use the following test to modify V so that all primes in 
V satisfy (Ic) and (Id) with high probability: 

deleteUnluckyPrimesNormal: Define an equivalence relation on V by setting 
pr^q :^=^ (LM(G(p)) = LM(G(g)) and LM(i?(p)) = LM{R{q))). Then replace V 
by an equivalence class of larges^ cardinality, and change QV and TZV accordingly. 

Only now, we lift the Grobner bases in QV and TZV to sets of polynomials G and 
R, respectively. Since we do not know whether all primes in the chosen equivalence 
class are indeed lucky and whether the class is sufficiently large, a final verification 
step is needed: We have to check whether ^ (G) is integral over A and normal. Since 
this can be expensive, especially if the result is false, we test the result at first in 
positive characteristic: 



If applicable, take Remark 4.7 below into account 
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pTestNormal: Randomly choose a prime number p ^ V such that Ap is a do- 
main, dp is a non-zero element in the conductor of Ap, and p does not divide the 
numerator and denominator of any coefficient occurring in a polynomial in G, R, 
or {/i, . . . , fr}- Return true if -j- (Gp) is the normalization of Ap and satisfies the 
relations Rp, and false otherwise. 

If pTestNormal returns false, then V is not sufficiently large for A or not all 
primes in V are lucky (or the extra prime chosen in pTestNormal is unlucky). 
In this case, we enlarge the set V hy t primes not used so far and repeat the 
whole process. On the other hand, if pTestNormal returns true, then most likely 
G = G{Q) and, thus, ^ (G) = A. It makes, then, sense to verify the result over the 
rationals by applying the following lemma. If the verification fails, we enlarge V 
and repeat the process. 

Lemma 4.6. With notation as above, the ring ^ (G) C Q{A) is the normalization 
of A if and only if the following two conditions hold: 

(1) The ring ^ (G) is integral over A. This holds if G and R are Grobner bases, 
and the elements of 2G satisfy the relations R. 

(2) The ring ^ (G) is normal. Equivalently, ^ (G) satisfies the conditions of 
the Grauert and Remmert criterion. 

Proof. If 2 (G) is integral over A, then ^ (G) C A. If ^ (G) is also normal, then 
equality holds. Note that if i? is a Grobner basis, then dim {R) = dim {R{p)) for 
all p & V. Hence, if the elements of gG satisfy the relations R, and G is a Grobner 
basis, then ^ (G) is integral over A. □ 

We summarize modular normalization in Algorithm [3j 

Remark 4.7. If the loop in Algorithm [3] requires more than one round, we have to 
apply deleteUnluckyPrimesNormal in Step 12 with some care. Otherwise, it 
may happen that always classes containing only unlucky primes are selected. To 
avoid this problem, when determining the cardinality of the classes considered in a 
certain round of the loop, we count all prime numbers in the class selected in the 
previous round as just one element. Then V will eventually contain lucky primes 
and termination of the algorithm is ensured by Lemma |4.5| and Remark |4.2| 

Remark 4.8. In Algorithm jsj the normalizations {G{p)) can be computed in 
parallel. Furthermore, we can parallelize the final verifications of integrality and 
normality. 

Remark 4.9. Algorithm [s] is also applicable without the final tests (that is, without 
the verification that ^ (G) C Q(^) is integral over A and normal). In this case, the 
algorithm is probabilistic, that is, the output ^ (G) C Q{A) is the normalization of 
A only with high probability. This usually accelerates the algorithm considerably. 

Remark 4.10. The computation of the algebra structure R of the normalization 
via lifting of the relations R{p) may require a large number of primes. Hence, if 
the number of cores available is limited, a better choice is to obtain just G by the 
modular approach and then compute the relations i?(0) over the rationals. For this 
approach, the initial ideals of the relations need not be tested in deleteUnluck- 
yPrimesNormal and pTestNormal. 
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Algorithm 3 Modular normalization 
Input: A prime ideal I CQ[X]. 

Output: A Grobner basis G C Q[X] and d e Q[X] such that ^ (G) C q{A) is the 
normalization oi A = Q[X]/I. 

1: compute M, the Jacobian ideal of A; 

2: choose a polynomial d e Q[X] representing a non-zero element d S M; 
3: choose V, a list of random primes; 

4: gv = (d,nr = 0; 

5: loop 

6: for p eV do 

7: if Ap is not a domain or dp G Ap is zero or dp is not contained in the 

conductor of Ap then 
8: delete p; 

9: else 

10: use the GLS algorithm to compute G{p), the reduced Grobner basis 

such that -J- (G(p)) C Q{Ap) is the normalization of Ap, and R{p), the 
reduced Grobner basis of the ideal of algebra relations; 

11: gv = gvu {Gip)}, nv^nvu {Rip)}-, 

12: {gV,nV,V) = DELETEUNLUCKYPRIMESNORMAL(^;'P,7^'P,7'); 

13: lift {gV,nV,V) to G C Q[A] and R C W[Ti,. . . ,T,] via Chinese remain- 
dering and the Farcy rational map; 

14: if the hft succeeds and pTestNormal(J, d, G, R, V) then 

15: if \ (G) C Q(A) is integral over A then 

16: if \ (G) C Q(A) is normal then 

17: return {G,d); 

18: enlarge V; 



5. Timings 

We compare the GLS normalization algorithirj^ (denoted in the tables below 
by normal) with Algorithm [2] from Section |3] (locNormal) and Algorithm [s] from 
Section [4] (modNormal|^ For all modular computations, we use the simplified algo- 
rithm as specified in Remark |4.10[ Note that at this writing, modularized versions 
of locNormal have not yet been implemented. 

In many cases, it turns out that the final verification is a time consuming step 
of modNormal. To quantify the improvement of computation times by omitting 
the verification, we give timings for the resulting, now probabilistic, version of 
Algorithm [3] (denoted by modNormal* in the tables). In all examples computed so 
far, the result of the probabilistic algorithm is indeed correct. 

All timings are in seconds on an AMD Opteron 6174 machine with 48 cores, 
2.2 GHz, and 128 GB of RAM, running a Linux operating system. Computations 
which did not finish within 2000 seconds are marked by a dash. The maximum 
number of cores used is written in square brackets. For the single core version of 
modNormal, we indicate the number of primes used by the algorithm in brackets. 



We use the implementation available in the Singular library normal . lib. 
^To implement our algorithms, we have created the Singular libraries modnormal . lib and 
locnormal . lib. 



PARALLEL ALGORITHMS FOR NORMALIZATION 



16 



So far, in Singular, the computation of associated primes via fast modular 
methods has only been implemented for the zero-dimensional case. As the compu- 
tation of associated primes is required by the local approach, we first focus on the 
case of curves, where the singular locus is zero-dimensional. 

The projective plane curves defined by the equations 

fik = {X''+^ + Y''+^ + Z^+^)^ - 4 l^x^+^Y^+'^ + Y^+^Z^+^ + z''+'^X''+^) 
were constructed in Hirano (1992)) . They have 3 (A: + 1) singularities of type A^, 



provided that k is even. If k is odd, the curves are reducible, in which case the 
normalization algorithms still work in the same way as in the irreducible case as 
long as they do not detect a zerodivisor. After the coordinate transformation 
Z I— 3^ — 2Y + Z, all singularities of the projective curves lie in the afhne chart 
Z 0. We apply the algorithm to the affine curves. The timings for k = 2, ...,5 
are shown in Table [T] 

Table 1. Timings for plane curves with many singularities: 





fl,2 


/l,3 


flA 


fl,5 


normal [1] 


.34 


14 






locNormal [1] 


.57 


2.0 


2.1 


38 


locNormal [20] 


.42 


1.3 


1.4 


11 


modNormal [1] 


4.4(3) 


73 (4) 


250(5) 




modNormal [10] 


4.1 


68 


240 




modNormal* [1] 


.57(3) 


7.4(4) 


11 (5) 




modNormal* [10] 


.31 


2.1 


2.5 





Both the local and the probabilistic modular approach have a better performance 
than the GLS algorithm, and they improve further in their parallel versions. The 
modular algorithm with final verification is slower, but can still handle much bigger 
examples than GLS. 

Timings for the affine plane curves defined by 

h,k = ({X - 1)'' - y'^)({x + 1)*-' - y'')(x'' - Y^){{x - 2)*^ - y^)((x + 2)*^ - y^) + y^'\ 

/g = X"' + y"' + (X - 2y + 1)^" + 2{X^'{X -2Y + if - X-'Y-' + Y-'{X -2Y + if), 

/4 = (y5 _,_ 2x*)(y^ + 7{x - i)'')({y + 5)^ + 2X^2) + y", 

/s = 9127158539954X"' + 3212722859346X*y^ + 228715574724X''y'' 

- 34263110700X*'y'^ - 5431439286X2y® - 201803238y^'' - 134266087241X* 

- 15052058268X'^y2 + 12024807786X''y* + 506101284X2y^ - 20217284iy* 
+ 761328152X'^ - 128361096X''y2 + 47970216X^y* - 6697080y'' 

- 2042158X* + 660492X^y2 - 84366y'' + 2494X^ - 474y2 - 1, 
are presented in Table [2] 

In Table H we consider surfaces in cut out by 
fgj. = XY(X - YfX + Y)(Y - 1)Z + (X*" - y2)(X^° - {Y - if), 
fyj. = - (y2 - 1234X^)'= (15791X2 - Y^){12ZIY'^ - X'\X + 158))(1357y^ - 3X"), 
/g = - ((13X - 17y)(5X2 - 7y^)(3X^ - 2y2)(19y2 - 2ZX'^{X + 29)))^ 
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Table 2. Timings for plane curves with various types of singularities: 





/2,7 


/2,8 


/2,9 


/a 


/4 


/5 


normal [1] 


7.7 


12 


383 




474 


1620 


locNormal [1] 


4.4 


13 


118 


1.9 


19 


1.2 


locNormal[20] 


1.4 


3.3 


31 


1.4 


18 


.93 


modNormal [1] 


38 (3) 


69 (3) 


146(3) 


142(3) 




50 (8) 


modNormal [10] 


38 


69 


146 


84 




43 


modNormal* [1] 


.70(3) 


1.2(3) 


1.2(3) 


88 (3) 


9.8(3) 


7.0 (8) 


modNormal* [10] 


.47 


.70 


.74 


30 


4.7 


.98 



We omit the verification step in the modular algorithm, as this is too time 
consuming. 

Table 3. Timings for the normalization of surfaces in A"^: 





/e.ii 


/6,12 


fe,i3 


h,2 


/7,3 


/8 


normal [1] 


2.6 


11 


6.4 








locNormal [1] 


.25 


.26 


.29 


80 


113 


70 


locNormal [20] 


.21 


.22 


.24 


80 


113 


70 


modNormal* [1] 


2.2 (2) 


.60 (2) 


.78 (2) 


12 (5) 


17 (5) 


2.3(2) 


modNormal* [10] 


1.5 


.52 


.67 


3.5 


4.7 


1.7 



Note that the performance of the local approach will be considerably improved 
as soon as modular primary decomposition in higher dimension will be available in 
Singular. 

Timings for the curves in A'^ defined by the ideals 

/g = {Z^ - (19y2 _ 2-iX'^{X + 29))^ - (liy2 _ 13Z2(Z + I))'') 
and the surface in A''^ defined by 

/lo = - {Y^ - 123456V1^2)(15791X2 - Y^f, 
WZ - (123ir^ - X{lllX + 158))) 

are given in Table |4j 

Table 4. Timings for curves in A'^ and a surface in A^: 





^9,1 


l9,2 


ho 


normal [1] 


3.2 




150 


locNormal [1] 


4.2 


36 


83 


locNormal [20] 


4.1 


35 


82 


modNormal [1] 






28 (4) 


modNormal [10] 






14 


modNormal* [1] 


8.9(5) 




8.4(4) 


modNormal* [10] 


2.1 




2.5 



To summarize, both the local and the probabilistic modular approach provide a 
significant improvement over the GLS algorithm in computation times and size of 
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the examples covered. The probabiUstic method is very stable in the sense that it 
produces the correct result in all examples computed so far. As usual, the verifica- 
tion step in the modular setup is the most time consuming task, and a refinement 
of this step will be the focus of further research. The modular technique parallelizes 
completely, the local approach parallelizes best if the complexity distributes evenly 
over the minimal strata of the singular locus. In general, the localization technique, 
even when not run in parallel, is a major improvement to the GLS algorithm. Note, 
that the local contribution can also be obtained by other means. See, for example, 
Bohm et al. (2011)] for a fast method in the case of curves, using Hensel lifting 
and Puiseux series. 
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